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Abstract 

> 

Measurements of the shapes of scintillation pulses produced by nuclear recoils, alpha 

particles and photons in Nal(Tl) crystals at visible energies of 10-100 keV have been 



performed in order to investigate possible sources of background in Nal(Tl) dark matter 



experiments and, in particular, the possible origin of the anomalous fast time constant 
events observed in the UK Dark Matter Collaboration experiments at Boulby mine [lj. 
Pulses initiated by X-rays (via photoelectric effect close to the surface of the crystal) were 
found not to differ from those produced by high-energy photons (via Compton electrons 
Ph! inside the crystal) within experimental errors. However, pulses induced by alpha particles 

(degraded from an external MeV source) were found to be ~ 10% faster than those of 
nuclear recoils, but insufficiently fast to account for the anomalous events. 
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1. Introduction 

Several Nal-based experiments searching for Weakly Interacting Massive Particles 
(WIMPs), as possible constituents of galactic dark matter, use pulse-shape discrimina- 
tion to separate the nuclear recoil signals expected from WIMP elastic scattering from 
electron recoils due to gamma background (see for example |2|, |||, |J). The feasibility 
of this technique arises because pulses initiated by nuclear recoils in Nal(Tl) are known 
to be typically 30% faster than those due to electron recoils |3], [|, |J. In each case the 
integrated scintillation pulses can be adequately fitted by assuming an exponential decay 
of time constant r. The resulting distribution of r, can be approximated by a gaussian 
in ln(r). The shape of the pulses can then be characterised by the mean value of r, r , 
in the gaussian ||. The difference between r values for nuclear and electron recoils is 
of the order of 30% for measured energies above 20 keV, decreasing almost to zero at 2 
keV. This difference, divided by the distribution width, is a measure of the discrimination 
which can be improved by optimising the crystal operating temperature fjj. 

Assuming operation of the technique at energies higher than ~4 keV |l], || then, apart 
from gammas, other potential sources of background are neutrons and alpha particles. 
Neutron interactions are expected to yield events indistinguishable from WIMP inter- 
actions (hence neutrons are usually used to determine the response of Nal to nuclear 
recoils since neutron scattering on nuclei is similar to WIMP scattering in terms of nu- 
clear recoil generation). However, neutron background can be sufficiently suppressed by 
shielding the detectors from neighbouring radioactivity and using deep underground sites 
to avoid neutrons from cosmic-ray muons. In the case of alpha particles any background 
at low energies cannot arise from the likely small contamination of uranium and thorium 
in Nal(Tl) because the alphas from uranium and thorium decay chains have energies 
exceeding 1 MeV. However, high energy alphas from activity in surrounding materials 
potentially may penetrate into the crystal and deposit energy in the surface layers. Their 
path in the sensitive surface layer of a crystal, and hence their energy deposition, can then 
be small depending on the point of their production. Thus such alphas could be respon- 
sible for background events and may be responsible for the anomalous fast time constant 
events (pulses faster than recoil-like pulses) observed in UKDMC Nal(Tl) experiments 
undertaken at Boulby mine JT], [?]]. 

In this paper we report characterisation of the shape of pulses induced by neutrons, 
alphas and electron recoils in Nal(Tl) in order to investigate possible sources of back- 
ground in Nal(Tl) dark matter experiments. Similar measurements have been performed 
previously in Nal(Tl) (see, for example, [Q, 0, [5|] and references therein). However, pre- 
sented here are results obtained with a single Nal(Tl) crystal at low energies (10-100 keV), 
allowing direct comparison of the pulses induced by the different particles in the energy 
range of prime interest for dark matter experiments. 

2. Experimental set-up and analysis technique 

Measurements were performed using a Hilger Analytical Ltd. unencapsulated Nal(Tl) 
crystal of dimensions 2 inch x 2 inch diameter cooled, using a purpose-built copper cryo- 



stat, to 11°C with control so that variation did not exceed 0.5°C Further details of the 
apparatus can be found in || ||]. The temperature chosen was close to that of the crystals 
used in the UK Dark Matter Collaboration (UKDMC) experiments at Boulby mine |], |7j • 
The crystal was viewed through a pair of silica light guides by two 3-inch ETL type 9265 
photomultiplier tubes (PMTs). Integrated pulses from the PMTs were digitised using a 
LeCroy 9430 oscilloscope driven by a Macintosh computer running Labview-based data 
acquisition software identical to that used with the UKDMC dark matter detectors. The 
digitised pulse shapes were passed to the computer and stored on disk. For the final 
analysis the sum of the pulses from two PMTs was used. The apparatus was found to 
yield a total light output corresponded to ~5.5 photoelectrons/keV. 

Energy calibrations for the tests were performed with a 57 Co gamma source (122 keV). 
A 60 Co gamma source was used to obtain pulse shapes from high-energy gammas. These 
gammas undergo Compton scattering inside the crystal producing electrons which deposit 
energy measured by the detector. To measure the shapes of the integrated pulses induced 
by nuclear recoils (sodium and iodine), the crystal was irradiated by neutrons from a 252 Cf 
source. Neutron energies were decreased by shielding the source with a 4 cm thick lead 
block. All aforementioned sources were placed outside the copper vessel containing the 
crystal, light guides and PMTs. An 241 Am source was used to irradiate the crystal with 
alpha particles. However, since the path length of alphas does not exceed tens of microns 
the source was attached to the bare cylindrical face of the crystal inside the cryostat. To 
decrease the energy of the alphas down to the keV range of interest a few, thin (~10 /xm), 
layers of plastic were placed between the source and the crystal. The 60 keV gamma-line 
from the 241 Am source also allowed independent energy calibration of the detector and 
measurement of the shape of the pulses initiated by electrons near the crystal surface. 
Such electrons are produced via the photoelectric effect. 

Analysis was performed, using the procedure described in || , by fitting an exponential 
to each integrated pulse to obtain the index of the exponent, r - faster pulses (due to 
nuclear recoils, for example) having smaller values of r, than slower pulses such as from 
gammas. For each experiment histograms were generated of the number of detected 
events versus the value of r - referred to here as r-distributions. The r-distributions are 
known to be dependent on the measured energy. To reduce such dependence the energy 
range (0-100 keV) was subdivided into energy bins of 10 keV width. We note that the 
energy threshold was about 10 keV. As shown in || (see also references therein) the r- 
distribution for each population of pulses can be approximated by a gaussian in ln(r) (for 
a more detailed discussion of the distributions see || and references therein): 

dN N r -(lnr-lnr ) 

exp 
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The r-distributions were fitted with a gaussian in ln(r) with several free parameters 
as follows. In the case of events from the 60 Co gamma source a 3-parameter fit was used 
with free parameters r , w and iV . In the experiments with the 252 Cf neutron source 
both neutrons and gammas (from the source as well as from local radioactivity) were 
detected. The resulting r- distribution can thus be fitted with two gaussians. However, the 
parameters r 07 and w for events initiated by gammas (effectively by Compton electrons) 
are known from experiments with the 60 Co gamma source. Assuming the value of w (called 



the width parameter) for the neutron distribution (where the pulses are due to nuclear 
recoils) is the same as that of the gamma distribution, since the width is determined 
mainly by the number of photoelectrons, again a 3-parameter fit can be applied. In this 
case the free parameters are: number of neutrons, N on , number of gammas, iV 07 , and 
the mean value of the exponent for the neutron distribution, r on . In practice, for direct 
comparison of the gamma distributions obtained with different sources, the r-distribution 
of gamma events measured with the gamma source was used, instead of the gaussian fit, 
to approximate the distribution of gamma events measured with the neutron source. 

In experiments with the alpha source both alphas and gammas were detected. To 
approximate the resulting distribution, we again used the r-distribution of the gamma 
events measured with 60 Co gamma source and a gaussian fit to the alpha distribution 
with 3 free parameters: number of alphas, N oa , number of gammas, iV 07 , and the mean 
value of the exponent of the alpha distribution, r oa . Furthermore, by making use of the 60 
keV X-rays from the 241 Am alpha source it was also possible to evaluate t x, attributed 
to pulses due to X-ray events initiated via the photoelectric effect near the surface of the 
crystal. Finally, to compare the 4 populations of events (initiated by gammas, neutrons, 
alphas and X-rays), we compared the values of r 07 , r on , r oa and r x- 

3. Results and discussion 

Measured r-distributions for events in two example energy bins are plotted in Fig- 
ure la (30-40 keV, alpha source), lb (55-65 keV, alpha source), lc (30-40 keV, neutron 
source) and Id (55-65 keV, neutron source). Plus signs show the data collected with the 
aforementioned sources. Open squares correspond to the data collected with the 60 Co 
gamma source normalised using the best fit procedure. Dotted curves show the fits to the 
neutron (alpha) distributions. 

Data collected with the gamma source (squares in Figure 1) match well the right-hand 
parts of distributions obtained with the neutron or alpha sources (these parts correspond 
to gamma events detected in the experiments with the neutron or alpha sources). This 
is true also for the 55-65 keV range with the alpha source where the gamma (right-hand) 
part of the r-distribution is dominated by X-rays from the 241 Am source. The amplitude 
of the right-hand peak at 55-65 keV (Figure lb) is several times more than that at 30-40 
keV (Figure la) (note the logarithmic scale of the y-axes), showing the presence of the 
strong 60 keV line superimposed on the background due to Compton electrons. This 
means that the r-distribution, and hence the basic shape of the pulses due to the X-rays, 
does not differ from that of the Compton electrons initiated by high-energy gammas. It is 
clear also that the positions of the left-hand peaks in the experiment with the alpha source 
(Figures la and lb) are shifted to the left with respect to the positions of the left-hand 
peaks in the experiment with the neutron source (Figures lc and Id, respectively). This 
is an indication that r oa is less than r on and, hence, the pulses due to alphas are faster 
than the pulses due to nuclear recoils (note that it was shown in |5], |[ that pulses due to 
sodium recoils are indistinguishable from those initiated by iodine recoils at all energies 
of interest). 

Figure 2 shows the fits to the r-distributions from the gamma- (solid curve), neutron- 



(dashed curve) and alpha- (dotted curve) induced events for the energy bin 30-40 keV. 
The total number of events in each case is normalised to unity. 

The results are summarised quantitatively in terms of r in Table 1. The typical error 
in t is of the order of 2-5 ns (arising from the statistics of the fit), except for the first 
energy bin in the experiment with the alpha source where the error of T oa is 15 ns - being 
higher due to the smaller number of detected alphas. The value of t x, obtained from the 
fit to the right-hand peak of the r-distribution at 55-65 keV with the alpha source (see 
Figure 1), is 322±2 ns, in good agreement with the value of r 07 (2nd column of Table 1). 



However, this does not agree with the conclusion of ||10| , where the shapes of the pulses 
due to X-rays from an 241 Am source and Compton electrons from high-energy gammas 
were found to be different (the shape of the pulses due to X-rays was found to be similar 
to that of nuclear recoils). The values of r oa (4th column of Table 1) are on average 10% 
smaller than those of r on (3rd column of Table 1). 

The values of r are known to vary from one crystal to another depending on the 
growth technology, Tl doping, temperature and other factors ]I|, 0, [J[ [I], [5||. However, 
the ratios, for instance of r on to r 07 are known to be quasi-independent of the crystal 
for fixed energy and temperature (we found it to decrease from 0.80 down to 0.76 with 
increasing energy from 10 to 80 keV in the crystal with anomalous events, currently under 
operation in the Boulby mine). The ratios Ton/r^, r oa /r ^ and r oa /r on are shown in the 
5th, 6th and 7th columns of Table 1, respectively. The first two slightly decrease with 
increasing energy, while T oa /r on remains almost constant. The average ratio is < T oa /r on > 
= 0.90 ± 0.01. The ratio < r oa /r on > is higher than the ratio < T oa /r on > = 0.79 ± 0.04 
found for the anomalous fast events (pulses faster than recoil-like pulses) observed in the 
UKDMC experiment at Boulby mine [1], 0]. This suggest that the anomalous events are 
not produced by external high energy alphas degraded in energy by a non-scintillating 
layer of material, assuming the ratio < r oa /r on > does not depend on the crystal. 

4. Conclusions 

The form of the pulses initiated by gammas, alphas, nuclear recoils and X-rays have 
been analysed in terms of the mean value, r D , of the gaussian distribution of exponent 
indices (see eq. (1)). The value of r oX of events initiated by X-rays (using the 60 keV line 
from an 241 Am source) was found to be the same as that of events induced by Compton 
electrons from high-energy gammas. The values of r oa for alpha events are smaller (by 
~ 10% on average) than those of r on for nuclear recoils induced by neutrons. However, the 
ratio of r oa JT on is higher than the corresponding ratio for the anomalous events to nuclear 
recoil events observed in the UKDMC experiment. This suggests that the anomalous 
events are not produced by external high energy alphas degraded in energy by a non- 
scintillating layer of material, assuming the ratio r oa /r on does not depend on the crystal. 
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Table 1: Mean values, r , of gaussian distributions (see eq. (1)) for gamma- (column 2), 
neutron- (column 3) and alpha- (column 4) induced events, and ratios of the mean values 
(columns 5, 6 and 7). The typical errors are 2-5 ns (in r , see text for details) and 0.02 
(in the ratios). 
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Figure 1: r-distributions for events collected with different sources: a) alpha source, 
measured energy 30-40 keV; b) alpha source, measured energy 55-65 keV; c) neutron 
source, measured energy 30-40 keV; d) neutron source, measured energy 55-65 keV. Plus 
signs show the data collected with the aforementioned sources. Open squares correspond 
to the measured distribution of Compton electrons (from 60 Co events) normalised using 
the best fit procedure. Dotted curves are the best fits to the left-hand parts of the 
distributions (alpha or neutron induced events). 
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Figure 2: Fits to r-distributions for gamma- (solid curve), neutron- (dashed curve) and 
alpha- (dotted curve) induced events for the energy bin 30-40 keV. The total number of 
events for each fit is normalised to unity. 



